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Abstract 



The highly structured search algorithm proposed by Hogg[Phys. Rev. Lett. 



80,2473(1998)] is implemented experimentally for the 1-SAT problem in a 

O 

r — ' single search step by using nuclear magnetic resonance technique with two- 



be readily extended to solving 1-SAT problem for more qubits in one step if 
the appropriate samples possessing more qubits are experimentally feasible. 
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Quantum computers[l-5] can outperform classical ones owing to utilization of quantum 
mechanical effects. The potential for quantum parallel computing mainly resulted from su- 
perpositions and entanglements of quantum states has made quantum computers capable 
of solving classically intractable problems (such as factoring large integers[3]) and finding 
tractable solutions more rapidly(e. g., searching an unsorted database[6]). However, car- 
rying out an actual quantum computing is much difficult although the rudiment of the 
quantum computing is easily understood and several quantum algorithm have been discov- 
ered. Among the approaches meeting the requirement for performing quantum computing, 
the nuclear magnetic resonance (NMR) technique [7, 8] is considered to be one of current 
promising method, with which the Deutsch-Jozsa algorithm[9] and Grover's search one [6] 
has been experimentally implemented with two-qubit samples [10- 13]. Recently, it was re- 
ported that an algorithmic benchmark has been experimentally realized using NMR with 
seven qubits[14]. It was easily found from those works that NMR can act as a considerably 
satisfactory small-scaled quantum computer to test the simple cases of the various quantum 
algorithm. In this contribution, also using NMR scheme with two-qubit sample, we report 
the first experimental demonstration of another quantum search algorithm. This algorithm, 
called a highly structured quantum search algorithm[15], was claimed to be much more ef- 
ficient than Grover's algorithm [6] that ignores the problem structure, and can be exploited 
in many practical examples instead of for academic study alone as an artificial problem[9]. 

The Hogg's algorithm is associated with the combinatorial searching known as nonde- 
terministic polynomial-time problems(NP)[16], an important class of intractable problems. 
Its prototype is the satisfiability problems(SAT)[16]. A SAT consists of a logic formula in n 
boolean variables, Vi,...,V n , and the requirement to find an assignment, specifying a value 
for each binary variable, that makes the formula true. The logic formula can be expressed 
as a conjunction of m clauses and each clauses is a disjunction of some variables. When 
all the clauses have exactly k variables, the problem is called k-SAT, with 1-SAT being the 
simplest case. Combinatorial search aims at seeking the solution of a SAT from the total 
of 2 n assignments. In general, the computational cost of solving a SAT grows exponentially 



with n in the worst case, making the SAT one of the most difficult NP problems[16]. For 
a few simple cases like 1-SAT and 2-SAT, however, solutions can be quickly and accurately 
determined without exponential cost using the regular structure of the problem. As each 
false clause for a given assignment is counted as a conflict, solutions are assignments with no 
conflicts. The classical search cost for such cases is O(n) due to each clause eliminating one 
value for a single variable, but the quantum counterpart with Hogg's algorithm for 1-SAT 
and maximally constrained k-SAT a single step [15], with the help of quantum parallelism 
and interference combined with the problem structure, which shows very high efficiency for 
large n. 

Actually, Hogg's algorithm consists of three stages[15]. First, one makes an equal 
superposition of bases, \s >, of the n two-state quantum system as the initial state 
\%pi >= 2~2 J2s I s >> with bit strings s being all 2 n assignments of n variables. Secondly, a 
transformation R is applied on \ifii >. The elements R ss of the diagonal matrix R are deter- 
mined by the conflicts c s of s, which are governed by the constraints expressed in clauses of 
the logic formula.lt is explicitly expressed as 

{\/2 cos (2c — 1)| for even m 
i c for odd m. (1) 

Finally, an operation U with the form 

{n— 1 
2 ~ cos [(n — m + 1 — 2d)j] for even m 
2 _! e i(n-m)f(_^ for odd m, 

is acted. It can be seen that the matrix elements U rs depend only on the Hamming distance 
d rjS between r and s. As it describes the inconsistency of all the corresponding bits of bit 
strings r and s, d rs has some correlation with conflicts c s . By exploiting the intrinsic relation 
between d r ^ s and c s and ingeniously constructing K ss and U rs with d r)S and c s , making the 
combined transformation UR on \ipi > once is able to find the expected state, i. e., UR 
\ipi >= \ipo >, which has equal amplitudes among solutions and no amplitudes among 
nonsolutions. 



In the view of quantum physics, Hogg's algorithm subtly makes use of the properties of 
quantum superposition and interference. When the calculating system is prepared on the 
equal superposition state, \ipi >= 2"? J2s I s >> a ^ the possible values, which can be assigned 
to the logical formula, regardless correct or incorrect, are embodied in the state. Then the 
state is subjected to the transform R, which gives a varied phase to each component state 
of the superposition state according to its conflicts c s . The final operation U destructively 
annihilates every component state that its conflict is not zero, construtively enhances the 
correct component state that its conflicts is zero, finds out the exact component state that 
satisfies the logical formula. In one word, the essence of searching achieved in a single step 
is to examine all the assignments simultaneously and to make use of quantum interference 
embedded in UR. By a single step we mean the number of search steps or sequentially 
examined assignments, not the number of elementary computational operations required. 

Our implementation of Hogg's algorithm was performed using solution NMR, applied 
to an ensemble of carbon-13 labelled chloroform molecules ( 13 CHCl3). Two heteronuclear 
spins X H and 13 C in 13 CHCl3 work as two qubits in our quantum computing. Similar to 
other experiments [9- 12, 14, 17-19] of quantum information processing using NMR, all the 
logic operations were realized by a specified sequence of radio-frequency (RF) pulses and 
the spin-spin coupling between the nuclei, and should be finished in the period of time 
much shorter than the relaxation times T x and T 2 of the nuclei in order to minimize the 
decoherence effect. 

A quantum circuit for achieving this highly structured search in a two-qubit system is 
shown in Fig. 1. The algorithm starts in the state \ip s >= |00 >, labelling the states of 
X H and 13 C spins from left to right respectively. As the state of a bulk sample at room 
temperature is described by a density matrix p t ^ for a thermally equilibrated system,we 
prepared an effective pure state p — |00 >< 00 1 by temporal aver aging [20], with which a 
deviation density matrix with diagonal elements of diag(p s )=[l,0,0,0] could be extracted 
from the summation over cyclic permutations of the populations of the |01 >, |10 >, and 
1 11 > states three times on the thermal matrix p t h- The Hadamard gates H in Fig. 1 were 



implemented by using RF pulse sequence ( |j (tt) x or its equivalents and thus a uniform 
superposition state denoted by the density matrix pi = |^>j >< i/ji\ with diag(/9j) = [l, 1,1,1] 
was obtained. Then two crucial operations R and U were applied, which are related to the 
problem structure and expressions for the logic formula and lead to the enhancement of 
solution states and the cancel of nonsolution states. 

For our two-qubit system(i. e., n — 2), m, the number of clauses in the 1-SAT formula 
with satisfiable solutions can be assumed 2 or 1. When m = 3 and 4, no satisfying assignment 
exists due to over-constrained formula being unsatisfiable. Four satisfiable formulas in the 
m = 2 case are V± A V 2 , V\ A V 2 , V± A V 2 andVi A V 2 , where Vj (i — 1, 2) are the negation 
of V,. The corresponding solutions are easily found to be |11 >, 1 10 >, 1 01 > and 1 00 > 
respectively, where the logic variable with subscript 1 corresponds to the high qubit and 2 
the low qubit. Similarly, the formulas for m—1 are V 2 , V±, V 2 and V±, corresponding to the 
unnormalized answer states |01 > +|11 >, 1 10 > +|11 >, |00 > +|10 > and |00 > +|01 >, 
respectively. According to the number of states in the solutions, they can be called as single- 
and multiple- item searching. The algorithm on two instances of 1-SAT was experimentally 
demonstrated: one for m = 2 with clauses V\ and V 2 , and the other for m—1 with clause 
V 2 . 

First of all, we have to derived the expressions for R and U from Ref. [15]. With Eq.(l), 
after calculating the conflicts c of each string assigned to the logical formula respectively, 
the elements of the diagonal matrices Ry lA v 2 an d Rv 2 were calculated to be diag(Rv 1 AV 2 ) = 
[— 1, 1,1,1] and diag(Rv 2 ) = [i, 1, i, 1]- U could be represented by VFrVF, where W = H\®H2 
with Hi being the Hadamard gates, and T is a diagonal matrix with the diagonal elements 

{a/2 cos \{m — 2|r| — 1)| for even m 
j\r\ e -im^ £ or Q ^ m ^ 

Elements of the diagonal matrices T for m = 2 and m = 1 were evaluated respectively as 
diag(T m=2 ) = [1,1,1,-1] and diag(F m=1 ) = [l,i,i,-l]. 

In order to realize the transformations with NMR, we designed the following RF pulse 



sequences respectively: 

R Vi av 2 - (f) w (f) xi (f)_ w - ^ - (f) w (f) x2 ^.^ 

r ™=2 ; (f) ?/i (I )_ X1 (i)_ wi - ir - Cf) OT (i).^ (i)_ w 

r ™=i : (2 J^ (fj^ (I J yi V2)y 2 \V_ X2 \V_ y2 

To make the best use of the available coherence time and to diminish errors due to 
the increased number of RF pulses, we finally optimized or reduced the pulse sequence 
constructed from the executing pulses of H, T and R mentioned above with the help of 
NMR principle[21]. The reduced pulse sequence for (UR)v lA y 2 and (UR)y 2 shown in Fig. 2 
were applied upon p iy resulting the output matrix p = \ip >< ip$\ = (U R) piiJJ R)^ . 

In order to read out the results of Hogg's algorithm accurately, a method of quantum state 
tomography[22] was adopted to reconstruct all the elements of the output density matrix po- 
For this end, we applied different read-out pulses immediately after the searching pulses UR, 
one for each run. Explicitly, nine pulses EiE 2 , Ei (|J , Ei (|J ,(f) E 2 , (f J (f J , 
(f) (f) '(f) E 2> (f) (f) and (f ) (f) were exploited, where Ei (i = 1,2) means 
no pulse acted on the i-th spin. Then we calculated the signal intensities by integrating 
the proton and carbon spectral lines acquired in each run and reproduced the whole matrix 
Po using the least-square fitting to the coupled equations connecting the signal intensities 
with the matrix elements. The recovered matrices po from the experimental data were 
depicted in Fig. 3, along with the theoretical ones for comparison. It is clearly seen from 
Fig. 3 that only the elements of < ll|p (V\ A V 2 )\U >, < 01|p (l / 2 )|01 >, < 01|po(V 2 )|H >, 
< 1 1 1 po ( ^2 ) 1 1 > and < 1 1 1 p ( ^2 ) I H > are significant, which justifies the searched items to 
be 1 11 > for Vi/\^2 and unnormalized |01 > +|11 > for V 2 formula. Due to experimental 
errors, however, several elements that should be zero theoretically still have small amounts 
of modulus, with the maxima being 8% for single-item and 19% for multiple-item searching. 

Errors in the experiments result from several sources. The effective pure state |00 > 
prepared by temporal averaging is generally not ideal, and small residual populations in the 
states 1 01 >, 1 10 > and |11 > of the deviation density matrix will surely cause errors in the 
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later stage of experiments. The inhomogeneity of RF fields and static magnetic fields and 
imperfections of the pulse-length calibration are also important sources of errors. Although 
we have tried to minimize these effects by careful adjustments of and measurements with the 
apparatus, some factors were still not well under control. Lastly but not insignificantly, losses 
of the coherence and populations of the density matrices in the whole process of experiments 
will lead to deviations of measured data from those expected by theory, especially for the off- 
diagonal elements. Perhaps that is why the maximum error in the multiple-item searching 
is larger than that in the single-item case. In all senses, minimizing experimental errors in 
all available ways is essential to accurately operating NMR quantum computers. 

We have experimentally demonstrated a quantum algorithm for the highly structured 
combinatorial searching and found solutions to the 1-SAT problem in one algorithmic step 
almost surely using an NMR quantum computer with 2 qubits. By contrasting with other 
search methods that ignore the problem structure requiring 0(2") steps classically and 0(2 2 ) 
steps on quantum computers [6], the searching with Hogg's algorithm can be accomplished 
in principle more efficiently However, there is no difference between Hogg's algorithm and 
Grover's one when only two qubits are considered. The potential of the high efficiency owned 
by Hogg's algorithm will be displayed in solving various 1-SAT problems for n > 2 qubits 
in one step with current NMR technique, which can be readily extended from the present 
experiment, if the appropriate samples possessing more qubits can be experimentally feasible. 
Furthermore, Hogg's algorithm for 1-SAT could be generalized to the maximally constrained 
k-SAT problems[15], which have practical examples such as scheduling, finding low energy 
states of spin glasses and proteins, and automatic theorem proving. While scaling up an 
NMR computer to much large systems poses daunting challenges, building such a device 
with a few qubits by some creative approaches for demonstrating algorithms stated above 
and the others is promising. 

We thank Xijia Miao for help in the early stage of experiments. This work was supported 
by the Chinese Academy of Sciences and National Natural Science Foundation of China. 
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Captions of the figures 

Figure 1 A quantum circuit for implementing a highly structure search algorithm 
on a two-qubit computer. Two Hadamard gates H®H transform an effective pure states 
\ip s >= 1 00 > into a uniform superposition state \tpi >, which is then transformed to the 
answer state \i/jq > after the action of gates R and U. For the definition of R and U, see text. 

Figure 2 The reduced NMR pulse sequences used to execute a) (UR)y lA y 2 and b) 
(UR)v 2 operations. Narrow and wide boxes correspond to | and n pulses respectively. X 
and Y denote the pulses along the x- and y-axis, X and Y, opposite to the x- and y-axis. 
The time period r is set equal to jj, where J is the size of the spin-spin coupling between 
nuclei 1 H and 13 C. Experiments were performed using a Bruker ARX500 spectrometer. 
13 C-labelled CHCI3 were obtained from Cambridge Isotope Laboratories. 

Figure 3 Experimentally recovered and theoretically expected deviation density matri- 
ces, po, after completion of the combinatorial search with the aim of satisfying logic formulas 
of both Vi A V2 and V2. The ordinate represents the modulus of matrix elements of po (not 
normalized). The numbers 0, 1, 2 and 3 in the horizontal plane denote the subscripts 1 00 >, 
1 01 >, 1 10 > and |11 > of the elements, respectively, a) Experiment for V\ A V2, b) Theory 
for V\ A V2, c) Experiment for V2 and d) Theory for V2. 
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